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PACS numbers: 13.25.Hw, 13.60.Rj, 14.20.Lq 

I. Introduction 

The large mass of the B meson allows a wide spectrum 
of baryonic decays, which have, in total, a branching frac- 
tion of (6.8±0.6) % [T] . This makes B decays a good place 
to study the mechanisms of baryon production. One ap- 
proach to investigate the baryonization process in B de- 
cays is to measure and compare their exclusive branching 
fractions and study the dynamic structure of the decay, 
i.e., the influence of resonant subchannels. 

In this paper, we present a study of the decay B~ — > 
£c + Ptt~tt~ [2]. This decay is a resonant subchannel of 
the five body final state B~ — > A+pir + n~iT~, which has 
the largest hitherto known branching fraction among all 
baryonic B decays and hence is a good starting point 
for further investigations. The analyzed decay can be 
compared with B~ — > E®p~Tt + ir~ and B° —> U^ + pw~ , 
which have similar quark content and phase space. 

Large differences between the branching fractions 
of B~ E+ + piT~iT~, B~ -> £®pTT + ir~, and 

B° — > £^ + pir~ could indicate a considerable im- 
pact of intermediate states on baryonic B decays. 



For example the decay B~ — s- E®pir + ir~ allows 
a number of resonant three-body decays (includ- 
ing N, A , and p° resonances) that cannot occur 
in B~ — > S^ + pTi^ n~ . The importance of reso- 
nant subchannels can be quantified, e.g., by the ratio of 
[B{B~ -> S+ + p7r-7r-) + B{B~ -> E^ptt+tt-)} /B{B~ -> 

A+pTT + Tt~Tt~). 

The CLEO collaboration measured B(B~ — > 
A+pn+ir-iT-) = (22.5 ± 3.5 ± 5.8) x 1CT 4 and B{B~ -> 
S°pn + T:') = (4.4 ± 1.7 ± 1.1) x lfr 4 0. The decay 
B° — s- £+ + pir~ was measured by the CLEO [Hj and 
the Belle |3] collaborations. The Particle Data Group 
has calculated an average of B(B° —> E+ + pTr~) = 
(2.2 ± 0.7 ± 0.6) x 10~ 4 [I]. For all these branching frac- 
tions the first uncertainty is the combined statistical and 
systematic error and the second one is due to the uncer- 
tainty in B( A+ -> pK-Tr+) = (5.0 ± 1.3) % PQ. 

II. The BABAR experiment 

This analysis is based on a dataset of about 426 fb -1 , 
corresponding to 467 x 10 6 BB pairs. The sample 
was collected with the BABAR detector at the PEP-II 
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asymmetric-energy e + e~ storage ring, which was oper- 
ated at a center-of-mass energy equal to the T(AS) mass. 
For generation of simulated data (Monte Carlo - MC ) 
we use EvtGen [B] for event generation and GEANT4 [7] for 
detector simulation. 

The BaBAR detector is described in detail elsewhere 
[S]- The selection of proton, kaon, and pion candidates 
is based on measurements of the energy loss in the sil- 
icon vertex tracker and the drift chamber, and of the 
Cherenkov radiation in the detector of internally reflected 
Cherenkov light [H] . The average efficiency for pion iden- 
tification is approximately 95%, with a typical misiden- 
tification rate of 10 % due to other charged particles such 
as myons and kaons, depending on the momentum and 
the polar angle of the particle. The efficiency for kaon 
identification is about 95 % with a misidentification rate 
less than 5 % due to protons and pions. The efficiency for 
proton and antiproton identification is about 90% with 
a misidentification rate around 2 % due to kaons. 
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III. Decay reconstruction 



The decay B~ 
subchannel 1 



-> £+ + pir 7r is reconstructed in the 
-t A+tt + , A+ -> pK^n + . For the 
reconstruction of the B candidate the entire decay tree 
is fitted simultaneously. A vertex fit is performed for B~ , 
E^ + , and A£, and the x 2 fit probability is required to 
exceed 0.1 %. 

In order to suppress background, the invariant mass 
of the pK~ir + combination is required to satisfy 
2275MeV/c 2 < m pK -^+ < 2296 MeV/c 2 , i.e., compati- 
ble with coming from the decay A£~— > p K~ ir + . This 
selection corresponds to 2.8 times the observed width 
of reconstructed A£ candidates which are centered at 
m P K-TT+ = 2285.4 MeV/c 2 . The separation of signal 
from background in the B-candidate sample is obtained 
using two kinematic variables, AE = E^ — \fs/2 and 

where 



™es 



v /( S /2 + Pi .p B ) 2 /i? i 2 -|pB| 



is is the 



center-of-mass (CM) energy of the e + e~ pair and E^ the 
energy of the B candidate in the CM system. (Bi,p;) is 
the four-momentum vector of the e + e~ CM system and 
Pb the B-candidate momentum vector, both measured 
in the laboratory frame. For correctly reconstructed B 
decays, toes is centered at the B meson mass and AE 
is centered at zero. Throughout this analysis, B candi- 
dates are required to have Toes within 8MeV (3. 4c) of 
the measured B mass of toes = 5279.1 MeV. 

Figure [l] shows the distribution of AM = m( A+n + ) — 
m(A+) in data for candidates that satisfy the criteria 
described above for toes and m p K- n + and for which AE 
is between —60 MeV and +40 MeV. We perform a binned 
minimum x 2 fit using a second order polynomial for the 
description of the background and the sum of a Voigt 
distribution (the convolution of a Breit Wigner with a 
Gaussian) and a Gaussian to parameterize the E+ + sig- 
nal. A detailed explanation of the fit function is given in 
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FIG. 1: Fitted AM distribution for B candidates in data. 
All candidates are required to satisfy the selection criteria on 
m E s, nipx-^+ an <l AE. 
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FIG. 2: Fitted AE distribution in data with selection criteria 
applied to m pK - 7T + , tties, and AM. The goodness of the fit is 
X 2 /dof = 36/40. The AE signal region is between -60 MeV 
and 40 MeV, and enclosed by the two sideband regions that 
each have a width of 50 MeV. 



Section IV. The fitted E++ signal yield is V = 1020±95. 

Figure [2] shows the AE distribution in data for can- 
didates that satisfy the criteria described above for Toes 
and to p x-7t+ an d f° r which AM is between 0.157 GeV/ c 2 
and 0.178 GeV/c 2 . The latter is a selection of can- 
didates with an efficiency of 92% in signal MC. 

In the binned minimum x 2 fit we use t he sum of two 
Gaussian functions for the signal and a linear function for 
the background. The second Gaussian accommodates B 
decays with missing energy due to final state radiation. 
Each Gaussian has a mean parameter (/x) and a standard 
deviation (cr). The joint normalization is described by 



TABLE I: The parameters for the double-Gaussian function 
describing the signal contribution in the fit to the AE distri- 
bution shown in Fig. [2] fx is the fraction of the signal in the 
narrower Gaussian. 



Parameter 



Fit result 



iV sig 840 ± 55 

fx (70 ± 23)% 

Hx (-2.7 ± 1.2) MeV 

cri (11 ±1.6) MeV 

ll 2 (-16 ±14) MeV 

cr 2 (20 ± 5.6) MeV 
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iV s i g and the fraction of the first Gaussian is j\. We pa- 
rameterize the background shape of the AE distribution 
as a first-order polynomial which provides a good descrip- 
tion of the AE distribution for candidates in the toes 
sideband in the range 5.20 GeV/ c 2 < to E s < 5.26 GeV/ c 2 . 
All parameters are permitted to vary during fitting. Ta- 
ble |T] presents the resulting signal parameters. The signal 
yield is 840 ± 55 events. 



IV. Signal extraction 

There are two sources of background that contribute 
to the signal in AE and AM. The first one are B 
decays that have the same final state, in particular 
B~ — > yl+]5 7r + 7r~7r - , and the other one are B decays 
that have a among its decay products, e.g., B° 

— > S^ + p 7r~ 7T°. To reject this background we make a 
binwise fit using AM as discriminating variable to cre- 
ate a background-subtracted AE distribution from which 
we extract the true signal yield in order to determine 
B(B~ — > £+ + pir~ir~). The binwise fitting procedure is 
described in the following. 

After applying the selection in m pK - 7T + and Toes 
(no selection in AM), we divide the AE range 
(-105, 105) MeV into 14 equal slices and fit the AM 
distribution in each slice separately in the range 
0.14 GeV/ c 2 < AM < 0.2GeV/c 2 . In the fits the S++ 
signal is represented by the sum of a Voigt function and 
a Gaussian function. The Voigt distribution has four 
parameters (N S i g , fx, T, a) and models the signal peak re- 
gion, where \i is the mean of the Voigt and represents 
the + mass, which is fixed to the value obtained from 
an inclusive analysis of — > A+ir + candidates in the 
data. The parameter T is the intrinsic width of the ±7+ + 
and is fixed to the 2010 Review of Particle Properties 
(RPP) value [1], and a describes the detector resolution 
in AM for the E+ + determined, independently for each 
AE slice, from the signal MC. The remaining parameter 
N* [g is the fitted £+ + signal yield in each of the AE 
bins. 



FIG. 3: AE distribution for B~ -> X"++p 7r 7r candidates 
in data. The points with error bars represent the number of 
T++ candidates iV* ig from a fit to AM = m( Atir + )~m( At). 
All signal parameters for the AE distribution, except iV s ig, are 
fixed to those shown in Table [I] 



There is a correlation between AM and AE that is 
very prominent due to the inaccurate momentum mea- 
surement of the slow 7T + from the decay. As a result 
the E~t + signal has tails in the AM distribution that are 
modeled by the Gaussian function whose parameters are 
determined, independently for each AE slice, from the 
signal MC. The background is represented by a second- 
order polynomial. This shape was determined from the 
sidebands |A£;| € (50, 300) MeV and, compared to the 
other polynomials, gives the best \ 2 fit probability. The 
fits in AM determine the background level and the num- 
ber of Sf + baryons. 

Figure [3] shows the signal yield as a function of 

AE. We fit this distribution with the same functions de- 
scribed in Section II and fix the signal parameters, except 
for iV s ig, to those determined there. The true signal yield 
is N sW = 787 ± 43 events. 



V. Efficiency 

The efficiency is calculated from the simulated events. 
These events were generated uniformly in four-body 
phase space (PS), but the actual decay distribution is, 
a priori, unknown. Therefore, when calculating the ef- 
ficiency, we weight the MC events so that we reproduce 
the distributions of the two-body invariant mass distri- 
butions for the decay products of the B-candidates in 
data. The resulting efficiency is checked by repeating the 
procedure using the three-body masses, and then again 
using the angles between the B daughters in the B rest 
frame. The different procedures give an average efficiency 
of (11.3 ± 0.2( syst ))%, which is used to determine the 
branching fraction. Out of the efficiencies from the dif- 
ferent procedures, we use the maximum deviation from 



G 



the average efficiency as systematic uncertainty. The sta- 
tistical uncertainty, due to the use of the data, is negligi- 
ble compared to the statistical uncertainty in the event 
yield. The efficiency calculated using unweighted events 
is 11.0%. 



VI. Systematic uncertainties 

We estimate the uncertainty on the signal extraction 
by repeating the fit to AE in Figure [3] separately for 
each parameter while permitting this parameter to vary, 
by using a second order polynomial for the background, 
and by fitting only the background and subtracting its 
integral from the histogram content in the signal range 
in order to obtain an alternative signal yield. We add the 
absolute values of the deviations in the event yields from 
all of these variations. The resulting relative uncertainty 
on the signal yield is 4.0 %. Other systematic errors come 
from track reconstruction efficiency (2.4%) efficiency 
(1.8%), and the number of produced BB pairs in the 
data sample (1.1%). The total relative uncertainty on 
the branching fraction is 5.1 %. 

VII. Branching fraction results 

Using the results from the signal extraction, efficiency 
determination, and estimation of systematic errors we 
find 



E+ + pTT-n-) ■ B( A+ -> pK-n + ) = 



N sie 



BB 



(1.49 ± 0.08 (stat) ± 0.08 (syst) ) x 10- 5 , and (1) 



free from A* c + and A resonances and hence can be 
described by a four-body PS distribution. The symbol 
7r~ refers to the ir~ that has the lower momentum in 
the e + e~ CM system. The other tt~ is denoted as irj . 



resonances for B 
< m(Z++ir-) < 



We see no indication of A and yl* 4 
candidates in the range 3.050 GeV/c 2 < m( i c ' 1 n, 
3.450 GeV/c 2 , where the normalization of the PS distri- 
bution is determined by fitting the sideband-subtracted 
data (Fig. [4]). From the ratio of the efficiency-corrected 
integrals of the distributions in Fig|4j we calculate an 
effective PS fraction: 



B(B~ 



B{B- 



£++Ptt-tt~)ps _ 389 11.3% 

> S++PTT-TT-) ~ 11.0% ' 816 



(3) 

This percentage will be used to normalize the PS pro- 
jection in the two-body and three-body invariant mass 
distributions in Figs. [5]|6j and [7] 
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e ■ N BS • B( At -> pK-Tr+) 



= (2.98 ± 0.16 (stat) ± 0.15 (syst) ± 0.77 Ac ) x 10~ 4 . (2) 

In equation [2] the last error is due to the uncertainty in 
B{A+ ^ pK-n+). 

VIII. Fraction of PS distributed decays 

In order to compare the two-body and three-body in- 
variant masses of the B decay products in data with the 
prediction from PS, we determine an effective PS frac- 
tion of the total branching ratio. To do this, we assume 
that the resonant substructures are due to the interme- 
diate states A* c + — > £+ + tt~ and A — > pir~ , and 
the remainder is distributed according to four-body PS. 
We investigate all 2-D planes that are spanned by the 
two-body invariant masses of the B decay products, e.g. 
m( JC+ + 7r~) against m(p7r^), to look for a range that is 



FIG. 4: The m( St +7r a') distribution in data (points with 
error bars) and for simulated four-body phase space decays 
(histogram). The distribution in data results from a sideband 
subtraction in AE according to the definition in Figure [2] 



IX. Resonant subchannels 

Figure [5] shows the invariant mass distribution of 
p7r~ = {p7r~ , J57T^ } (sum of the distributions of m(pTT~) 
and m(pir^)) after sideband subtraction in AE (see Fig. 
[2] for the definition of the sidebands) and efficiency cor- 
rection. The efficiency correction here and in the other 
invariant masses of the B daughters is determined from 
PS MC for the particular mass that is considered. The 
differences between data and PS in the range m(p7r~) € 
(1.2, 1.7) GeV/c 2 are compatible with the existence of the 
resonances A — (1232, 1600, 1620). 

Figure [6] shows the invariant mass of = 
{E+ + tt~, Z , + + 7T f T} after efficiency correction and side- 
band subtraction in AE. The large number of events 
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FIG. 5: The m(pn~) distribution in data and simulated four- 
body phase space decays. The shaded vertical ranges repre- 
sent a width of one V and are centered at the average mass 
of A (1232) and A (1620), respectively. The parameters 
are taken from the RPP [1]. The range of A (1600) is not 
drawn since its parameters have large uncertainties. 



at threshold are consistent with the decay B~ — » 
yl c (2595) + p7r^. There are no significant signals for other 
A* + resonances. 
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FIG. 6: The m(X'+ + 7r _ ) distribution in data after efficiency 
correction and Ai5-sideband subtraction. The solid line shows 
four-body phase space decays. The shaded vertical ranges 
represent a width of one F and are centered at the average 
mass of the respective A* + resonance. The parameters are 
taken from the RPP pQ. 

In the three-body invariant mass distribution 
m( E+ + ir~n~) (Fig. [7| we see unexplained structures 
at 3.25GeV/c 2 , 3.8GeV/c 2 , and 4.2GeV/c 2 . However, 
due to the limited number of signal candidates, it is not 
possible to analyze these enhancements in more detail. 

We find no indication of a threshold enhancement in 
the baryon-antibaryon mass distribution. 



FIG. 7: The m( £+ + ir~n~) distribution in data and simu- 
lated four-body phase space decays. The histogram in data 
includes efficiency correction and A_B-sideband subtraction 
according to the definition in Fig. [2] 



X. Conclusions 

We have measured the branching fraction of B~ — > 
Ec + pir~ Tt~ . This improves on the previous measure- 
ment by CLEO [5]. Additionally, we find large devia- 
tions between data and the prediction from four-body 
PS in the distributions of the invariant masses of the 
B daughters. This may indicate the importance of res- 
onant substructures in baryonic B decays. The com- 
bined branching fraction of the resonant subchannels 
B~ — > E^ + p7r~TT~ and B~ — > E^pir^n^ makes up 
about 30 % of the branching fraction of the five-body 
final state B(B~ — > A+pir + ir~Tc~). 
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